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ON THE THEORY OF OSCILLATING AIRFOILS OF FINITE SPAN 
IN SUBSONIC COMPRESSIBLE FLOW 1 

By Ebic Reissneb 


SUMMARY 

The problem of the oscillating lifting surface of finite span in 
subsonic compressible flow is reduced to an integral equation. 
The kernel of ike integral equation is approximated by a simpler 
expression, on the basis of the assumption of sufficiently large 
aspect ratio. With this approximation the double integral 
occurring in the formulation of the problem is reduced to two 
single integrals, one of which is taken over the chord and the 
other over the span of the lifting surface. On the basis of this 
reduction the three-dimensional problem appears separated into 
two two-dimensional problems, one of them being effectively the 
problem of two-dimensional flow and the other being the problem 
of spanwise-circulation distribution. Earlier results concern- 
ing the oscillating lifting surface of finite span in incom- 
pressible flow are contained in the present more general results. 

INTRODUCTION 

The present report is concerned with the problem of the 
oscillating airfoil of finite span, within the frame of the linear- 
ized lifting-surface theory. The aim of this study is the 
development of a theory which incorporates simultaneously 
the effects of three-dimensionality of the flow and of com- 
pressibility of the fluid. As an exact solution of this problem, 
even within the limitations of the linearized theory, presents 
very great difficulties, it is worth while to work toward an 
approximate theory which is valid, provided the aspect ratio 
of the lifting surface is not too small. 

The author has previously obtained results of this nature 
for the case of incompressible flow (references 1 and 2). In 
this earlier work the known results for the problem of two- 
dimensional incompressible flow were contained as a special 
case. The present work generalizes these results so as to 
take account of compressibility in the subsonic range. Thus, 
the results of this report consist of a system of equations 
which contain, as special cases, both the author's results for 
the wing of finite span in incompressible flow and the results 
of Possio’s theory of two- dim ensional compressible flow 
(reference 3). 

The scope of the present results may briefly be described 
as follows. The starting point of the work is an integral- 
equation formulation of the problem of the lifting surface of 
finite span. The integrals which occur are double integrals 

* Supersedes NACA TN 1953, “On the Theory oT Oscillating Airfoils of Finite Span in 
Subsonic Compressible Flow” by Eric Reissner, 1949. 


and the functions to be determined are functions of two in- 
dependent variables. The essential step of the present work 
is to replace the actual kernel of the integral equation by an 
approximate kernel in such a way that the double integrals 
are reduced to single integrals over the range of either one 
of the two independent variables. In this. way the problem 
is reduced to two problems which are to be solved separately. 
The first of these two problems is of the same nature as the 
Possio problem of two-dimensional compressible flow. The 
second of these problems is of the same nature as the problem 
of the Prandtl lifting-line theory for the wing of finite span 
in uniform motion. 

As in the theory of incompressible flow, this reduction of 
the double-integral problem to two single-integral problems 
depends crucially on the assumption of sufficiently large 
aspect ratio. While “sufficiently large” aspect ratios might 
be thought to be aspect ratios of about 3, definite statements 
of this nature must be based on experimental evidence, as 
long as no exact solutions exist for the three-dimensional 
problem of the oscillating lifting surface in compressible flow. 

It is perhaps worth while to state explicitly that the pres- 
ent problem is quite different from the corresponding problem 
for supersonic flow. 

It might also be added that there are reasons to believe 
that it is not satisfactory, even approximately, to superimpose 
aspect-ratio corrections for incompressible flow and com- 
pressibility corrections for two-dimensional flow in order to 
obtain corrections for the combined effect. . This latter point 
is one of the reasons for the present study. 

This work was conducted at the Massachusetts Institute 
of Technology under the sponsorship and with the financial 
assistance of the National Advisory Committee for Aero- 
nautics. 

SYMBOLS 


X,Y,Z 

U 


t 

H 


u, v,w 


Po 


Cartesian coordinates 
main-stream velocity in X-direction 
time 

defined by equation of lifting surface 
Z=H(X,Y,t ) 

components of velocity change caused by pres- 
ence of lifting surface 

region in X,T-plane occupied by projection of 
lifting surface . 

density of stream flowing with velocity U 
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density and pressure changes, respectively, asso- 
ciated with velocity changes u, v, and w 
velocity of sound in main stream (a?=dp/dp 0 ) 
coordinate of trailing edge' of R a 
potential of velocity changes u, v, and w 
circular frequency of oscillation 
Mach number of main stream (17/a) 
real part of 

a length to be identified with the semichord 
of R a at midspan 

reduced-frequency parameter (aibfU) 
dimensionless coordinates, defined by equation 
(19) 

function defined by equation (21) 
parameter defined as 

parameter defined as K = jyyjp ’ 

k 

parameter defined as v= i_fi p 

region in £,y-plane corresponding to region R a 
in A,F-pIane 

coordinate of trailing edge of R a * 
coordinate of leading edge of R a * 
region in z,y-plane consisting of the strip to the 
right of the trailing edge of R a * 
entire a;,y-plane except for regions R a * and R a * 
function defined by equation (36) 
function defined by equation (37) 
variables of integration in accordance with 
equation (40) 
defined by equation (41) 

Hankel functions of second kind, and of zeroth 
and first order, respectively 
auxiliary variable of integration 
function defined by equation (51) 
function defined by equation (52) 
functions defined by equations (54), (55), and 
(56) 

order of magnitude of 
function defined by equation (38) 
function defined by equations (73) and (83) 
auxiliary variable of integration 
auxiliary variable of integration 
function defined by equation (76) 
ratio of semispan t o semic hord at midspan 
defined as 

local semichord divided by semichord b at mid- 
span 

a quantity indicating amount of sweep and 
defined by equation (86) 
dimensionless coordinate defined as 
x*=(x-x m *)/b* 

quantities defined by equation (87) 
quantities defined as 

g* OB*, y*) =9(x,y); x*(**, y*)=Kx, y) 

function defined by equation (89) 
function defined by equation (90) 


THE BOUNDARY-VALUE PROBLEM OF THE 
OSCILLATING LIFTING SURFACE 

It is assumed that a nearly plane, impenetrable surface is 
put into the path of an inviscid, flowing fluid which, except 
for the effect of this surface, possesses a uniform velocity V 
in the direction of the positive A’-axis. The impenetrable 
surface, henceforth called lifting surface, is taken to lie 
nearly in the A", F-plane, and its equation is written in the 
form Z=H(X,Y,t). When H— 0 no disturbance is caused. 
When the lifting surface is not exactly plane and parallel to 
the direction of U the velocity components (£7,0,0) are changed 
into (V-\-u,v,w) where u, v, and w depend on the form of the 
function H and on the shape of the region R a which is the 
projection of the lifting surface onto the A’/T-plnnc. 

The disturbances caused by the presence of the lifting 
surface are assumed to be small, in the sense that the differ- 
ential equations and boundary conditions of the problem arc 
linearized with respect to the disturbance velocity compon- 
ents u, v, and w and with respect to the pressure and density 
changes p and p caused by the presence of the lifting surface, 

Under these conditions the differential equations of the 
problem are the following: 


bu , TT bu d_ f p\ 

bt + U bX~ bX\ Po ) 

*>v , T t bu 5 fv\ 
bt ±U bX~ bY\pJ 

bw, Tr ^ W - — f Y-\ 

bt +u bX~ bZ \pj 

dp „ dp fbu bv,bw\_ n 
^ +U bX +pa \bX + W + bz)~ 0 


p=a 2 p (5) 

The quantity p 0 in equations (1) to (4) is .the density in the 
fluid flowing without disturbance and the quantity a in 
equation (5) is the velocity of sound in the undisturbed fluid; 
that is, a i =dp/dpi > . 

The boundary condition of no relative normal flow at the 
lifting sui'face is satisfied, within the frame of the linearized 
theory, instead of on the lifting surface itself, on tbc pro- 
jection of this surface onto the A'F-plane, 

X, Y inside R a , w =4|-f£/^ (6) 

The form of condition (6) (which holds on both sides of the 
lifting surface) indicates that w is an even function of Z, 
From equation (3) it follows then that p is an odd function 
of Z and the condition that the prossuro disturbance p is 
continuous, except when passing across the lifting surface, 
means that for Z—0 and 

X, Y outside R a , p — 0 (7) 

On the basis of conditions (7) and (6) the problem may be 

considered to consist in the determination of u, v, w, and p in 
the half space Z^0 with the boundary Z=0. 
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In addition to conditions (6) and (7) the following further 
conditions are prescribed in order to obtain an unambiguous 
solution. At the trailing edge of the lifting surface, 


Equation (10) becomes 


p = — Po 


U bXj 


(17) 


X = X t (Y ), p is finite (8) 

Finally, it is postulated as a condition “at infinity” that 
energy is traveling outward without reflection, in a manner 
to be defined more precisely in what follows for the case of 
simple harmonic motion. 


VELOCITY-POTENTIAL FORMULATION OF THE PROBLEM 

The problem as stated in equations (1) to (8) may be 
solved by means of a velocity potential <£, in terms of which 



v- 



y 


( 9 ) 



Combination of equations (9), (1), (2), and (3) results in the 
following expression for the pressure change p: 

*’=-'-(n+ p H) (10) 

Combination of equations (10), (5), (9), and (4) results in the 
following differential equation for <f>: 


The following dimensionless parameters and variables are 
now introduced: 



*=irj 


X 


X 

'b 


(18) 


„=yi=ip£ 


y 


(19) 


The differential equation (16) then assumes the form 


V s ?— i 


2 kAP £>? , PM 2 _ n 
l-APbx^l-AP* 


(20) 


Equation (20) is reduced further, for the purpose of elimi- 
nating first-derivative terms, by the following substitution: 


Choose 


APk 

M l-AP 


( 21 ) 

( 22 ) 



** , 


bX 2 ^ 

dF 2 ' 

' dZ* a 2 \ 


i*(t +cr sr )’* =0 


Equation (11) is to be solved in the half space Z^O subject 
to the following boundary conditions at Z=0: 


,, „ . . , „ b<j> bH . TT bE 
X, 1 inside B a , ^ + U $x 


( 12 ) 


X, Y outside R«, ^+U |j=0 


(13) 


and obtain as the equation for 


where 


vv+*v=o 

TcAl 

K ~1 -AP 


(23) 

(24) 


The expression for p and the boundary conditions must 
now be expressed in terms of the new independent variables 
x, y, and z and in terms of the new dependent variable 
From equations (17) and (19), it follows first that 


X=X t (Y), ^+U ||is finite 


(14) 



(25) 


and subj'ect to the condition of no energy reflection at in- 
finity. 

In what follows attention is restricted to the case of simple 
harmonic motion by writing: 

<t>(X,Y, Z,f)=e iw£ ?(Z,F,Z) (15) 

with corresponding expressions for H and p? 

Equation (11) now assumes the form 

Ip+If+S-^^s )'^ 0 (i6 > 


> It is perhaps not entirely superfluous to indicate that this is meant in the sense that, 
corresponding to a surface equation there is a pressure distribution 


Combination of equations (25) and (21) gives 

p= 

where v is defined as 

k 

v 1 —AP 


( 26 ) 

(27) 


The boundary condition (12) becomes 

x, y inside ^=-jJL= e -^(ikH+~) (28) 

where R a * is the region in the x, y-plane corresponding to 
B a in the X, F-plane. 
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The boundary condition (13) becomes 

x, y outside R*, iv\(/+^=0 (29) 

The condition of finite pressure at the trailing edge is now, 

x—x T (y), is finite (30) 

Finally, the condition of no energy reflection at infinity is 
written in the form 

4 ,!a J(x,y ) z)e- UT (31) 

where r 2 —x 2j ry 2J rz 2 and where / tends to zero as s tends 
to infinity. 3 

The boundary condition (29) is made more specific in the 
following manner. From condition (29), 

i{x,y,Q)=c{y)e~ irX (32) 

For any line y= Constant which does not pass through R*, 
it can be concluded from the condition of undisturbed flow 
at x=— oo that c(y)= 0. The same can be said for that 
portion of any line jf=Constant which is situated in front 
of the leading edge of the airfoil region, The situation is 
different for portions of lines Constant which do pass 
through R* and which are to the rear of the trailing edge 
of R*. The region to the rear of the trailing edge and 
bounded by lines y = Constant which are tangent toi? a * is 
called the wake region and is designated by R a *. The 

exterior of the region R a * and R a * is called the remaining 
region and is designated by R r *. Then c(y )= 0 in R* and, 
in general, c(y)^ 0 in R a *. 

Equation (29) is thus seen to be equivalent to. the follow- 
ing two equations: 


x , y inside R,*, 

(33 a) 

x, y inside R a *, \ p—c(y)e~ ir * 

(33b) 

If now a function A is introduced defined by 


A(y) = 2t[x T (y),y,0] - ' 

■ (34) 

then equation (33b) can be written in the form 


x,y inside R w *, ||=— yA (y)e ir < x T-*» 

(36) 


In view of equation (33a), there may also be written . 

' ( 36 > 

where x L (y) is the coordinate of the leading edge of the 
airfoil region R a *. 

SUMMARY OF THE RESULTANT BOUNDARY- VALUE PROBLEM 

Before proceeding with the solution of the problem as 
reduced in the foregoing section of this report its final 
formulation is recapitulated as follows. 

* This insures that £«/ as z tends to infinity and therewith that waves are 

traveling away from the source of the disturbance. 


Determine the solution of the differential equation 


vVH-*V=0 

(23) 

in the half space 2^0 subject to the following conditions 
at 2 = 0: 

x, y inside /?„*, ^<=g{x,y) 

(28a) 

x, y inside /?„ *, ||=— - y 

(35) 

x, y inside R r *, ||=0 

(32b) 

x=x T (y), || is finite 

(29a) 

and subject to the following condition at infinity: 


2— » oo t tp=jc~ tKT 

(31a) 

where f=o(z) and where r i =x 2 -{-y i -\-z i . 

Various quantities occurring in these equations are defined 
as. follows: 

;•• • s(z ’ v) -7rm( tkB+ w) 

(37) 

A (y)= \(x,y)dx 

(36a)* 

X(*,y)=2 

(38) 


The parameters k, h, and v are defined by equations (24), 
(22), and (27), respectively. The region R* follows from 
the airfoil region R a by multiplication with a scale factor 1/6 
in the r-direction and by multiplication with a scale factor 
■y/l — AFfb in the ^-direction. The region R u * is the strip 
extending from the trailing edge x—x T to 2 =°=, and the 
region R* is the remainder of the x, y-planc. 

The solution of the boundary-value problem is to be used 
to calculate the pressure-change amplitude p a at the lifting 
surface in accordance with the relation 



which follows from equation (26) . 

The solution of the problem as summarized will be 
approached through its reduction to an integral equation 
; for the quantity X as defined by equation (38). 

AN INTEGRAL REPRESENTATION FOR THE VALUES OF O./'/d* 

In this section it is proposed to derive a formula for the 
values of for the purpose of setting up the basic 

integral equation of the problem under consideration. To 
begin, results are taken which in essence are known and these 
are transformed in a way designed to facilitate. the subsequent 
transition from the exact double-integral equation of the 
problem to the approximate integral equation containing 
single integrals only. 
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The first formula is a representation of the values of 
bp/dx in the interior of the half space s^>0 in terms of the 
values of bfi/bz on the boundary 3=0 of the half space, as 
follows: 

In equation (40) and in ail that follows, 

r a =(a:— £)*+(y— v)*+z 2 (41) 


and the quantity X is according to equation (38) given by 
2 £>’/'(£, 7 , 0)/2>£. It is noted that this representation of bfjbx 
insures that the conditions at infinity as expressed in equa- 
tions (31) are satisfied. 

From equation (40), it follows that 


ay 
a xbz 


tS.h ’ ” ) 




(42) 


If it is now observed that the quantity r -1 e _i,tr is a solution 
of the differential equation (23) for rp, then equation (42) 
may be written in the alternate form 


ay 

a^a z 4 


7 / J X( *’ (^ + ^ +k2 ) (V) (43) 


In equation (43), it is noted that a 2 /a7/ 2 =a 2 /a?? 2 and the term 
in question is integrated by parts with respect to ij. In 
addition to this, by making use of the obvious identity 

x(£> i7)=M?j y)+[M£) v)~ y)J 


equation (43) is written as follows: 


bi&2 _ 4» J J x(e ’ (a? + “’)( 7 )&'- 

P/Wfci) MS, ?)](£,+ 



(44) 


In the first integral on the right of equation (44) the inte- 
gration with respect to rj may be carried out explicitly. 
When the remaining two integrals are absent, there is there- 
by obtained the appropriate form of the integral relation (44) 
which would follow if two-dimensional flow had been assumed 
from the beginning. 

The following formula expresses the integral in question 
in terms of a Hankel function 


f’“ 2 +z* 

V(*— 9*+(y— ■ 


dr, =ixH 0 ® [*V(*-0*-!-**l 


(45) 


The next step is to reduce the second integral on. the, right 
of equation (44) to an integral involving ax/dij. To this 
end, there is written 


y)l ^ dr , » £ + J“ (46) 
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The two integrals on the right of equation (46) are integrated 
by parts and the constants of integration are chosen in such 
a way that the integrated portions vanish. After some 
elementary transformations, this leads to the following 
formula : 



[\(£, ij)— X(£, y)] 


e — i* V(i— i,) a +2 a 

V(*-©*+( y-nY+z* n 


axkmir ffavfr-awi* Jfc l ^ 

^ y-n [J— V(*- ©*+«*+ r 2 J 


(47) 


Equations (45) and (47) are introduced into equation (44) 
and the following relation is obtained: 


ay - 

bxbz 2r 


[*V(*-0 2 +s 2 ]}dH- 




4xJ J bi)( y—v W 2 


+*■)/: 


— Ilf — «t g—ik V(±-£) 2 -r-z*+f 2 


d£-\- 


V(ar-S a +3*+f i 

a r rrtxi* «'-»>‘ 1),w. (48) 


The final step now consists in integrating equation (48) 
with respect to x between the limits — and x. The con- 
dition of undisturbed Sow far in front of the lifting surface 
makes (bip[bz) x= ~ a —Q. There is then obtained the relation 


ff X _ ^Ho C2) Mix'-Q'+z^dx^d^ 

J_ f pU iy— q| £> e-nVcr-a^-H: 8 

4 xj J dr, { y — r, V(iT— 4)*+S , + f* 


ly — u 1 
y—v 


■£[/: 

/: 


-lir-il Z'-W+P-H 3 


dflda:' 


V(z'-s) 2 +* 2 +r 2 




'+ 


d%dr j 


(49) 


Equation (49) expresses the values of- b\pfbz for 2>0' 
in terms of the values of X(f,-y)=25f(£,y I 0)/df and of 
5X (£,??) /dij. Because of this fact the present form of this 
integral relation is particularly well-adapted to separate 
three-dimensional from two-dimensional effects. 

THE INTEGRAL EQUATION OF THE OSCILLATING LIFTING 

SURFACE 


The integral equation of the problem is obtained by sub- 
stituting the information contained in equations (28), (32), 
: and (35) in equation (49) and by then letting z tend to zero. 
In so doing the integrals which contain infinities in the inte- 
grand are to be interpreted as Cauchy principal values. This 
latter step has been discussed in detail for the case of incom- 
pressible flow in reference 1. ISTo additional difficulties in 
this regard are encountered for the. present problem of com- 
pressible subsonic flow and the explicit justification for this 
step is therefore omitted in the present report. 
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The integral equation of the oscillating lifting surface in 
subsonic compressible flow is thus of the following form : 


1 J'-t(v) 

g(x,y ) = — ^ Jx L ( V ) V'ffl-b “ 0; «]<* ?+ 


A(y)e iKC T(v) f e-^K^x— §;ic]dZ+ 

ZtC Jx T (v) 

±^^G[(x-Q,(y- v y,K}didy- 

/ J ^ (y-y); 4 dl-dy 

The kernels K and Q are of the following form : 


(50) 


( | S - *!)]+ « 2 (i*'~ SD dx' (51) 


r . {y-y] *) * ... 

y-y \bx 

C x T C- l»-*i ( ,-f*V(j:'-e i, +r 2 "I ) 


V(*'- 8 *+f’ 

:' Jt {)*+(y=-T) 2 


+ 


C x 3 r e ‘-*«^C r ''*-8*+(v =! -f) 2 "1 

J-„ 3tj L V ’©*’4- Cy — ii) a J dx ( ' 52 ’ > 

Equation (50) holds when a: and y are inside R a * and is to 
be solved for X, in terms of A and g, where A and g are defined 
by equations (36) and (37), respectively. 

For the case of two-dimensional flow, b\fby=0 and the 
last two integrals in equation (50) are absent. The remainder 
of the present work has as its object the derivation of a 
procedure to take account of these last two integrals in an 
approximate manner which permits the calculation of the 
effect of Three-dimensionality of the flow in a simpler way 
than by actually solving the complete equation (50). In the 
derivation of this procedure the integral over the airfoil 
region R a * and the integral over the wake region R a * are 
treated separately. 

REDUCTION OF THE DOUBLE INTEGRAL OVER THE 
AIRFOIL REGION R a * 

Write . - — . 


where 


R.* 


0 f-ir— li g— «W i*-fl z +f 2 


'bxj- a voe-eM-r* 


di 




- vc^im 2 j 


, r* r p-lr 
f* b j" e -i*V(i'-© 2 +(y-i } 2 "J f 

3 


(54) 

( 55 ) 

( 56 ) 


When dealing with lifting surfaces of sufficiently elongated 
form, that is, with surfaces of sufficiently high aspect ratio, 
\y—y \> )>!*—?! over tire major portion of the surface. When 
this latter inequality holds, the terms J* can, as will be seen, 
be written in such a way that a relatively simple dominant- 
term can be separated in each of them. These dominant 
terms will be used for the approximation to be developed. 
There remains then the question concerning the validity of 
the approximation over that portion of R a * where the ine- 
quality (a:— £|<!<!h/— y\ does not hold. This question is 
answered as follows. It is assumed that bXfby varies suf- 
ficiently slowly with y so that in this portion of the region 
R a * it is effectively constant. If this is the case, all that is 
necessary is to take account of the fact that both the kernel 
G and the approximation to G to be obtained are odd func- 
tions of y—y, so that in both cases the contribution to the 
value of the integral coming from this portion of R a * can be 
neglected. 

The aforementioned argument is also implicit in the 
earlier derivations for the corresponding problem for in- 
compressible flow (references 1 and 2). There appears to 
be no reason to believe that this particular argument should 
be less applicable to the problem of subsonic compressible 
flow than to the problem of incompressible flown 
Equation (54) is written in the form 

r — If— vl [" x — £ 1 

1== ~J- (*-*m r Lv(T=^^F +iK(a:_e J rf! ' 

and this implies the following ordor-of-magnitude relations: 
T r-ty-ii e -**in / i , . \ , 


Equation (55) is written in the form 


(57) 


'■-jr 



-*jr 

[?«•««+/. 

dt (58) 

In equation (58), note that, when xis not too small and 
when \y— i?|»ja;— $[, 


/ e -i*Irl\ 

xwt-oijf) 

(59) 

r*-s 

Jo 

j-u^(z'> 2 -h 3 r "1 • 

=0 e T r (x J) 

Vmt* L irf 'J 

(60) 


Neglect the term in equation (60) compared with the term in 
equation (59) ; then, 




/: 


I* Mix 




( 01 ) 


Note, that this approximation ceases to be correct as k be- 
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comes smaller and smaller. However, when this is the case 
the contribution to the total coming from I 2 becomes negli- 
gible because of the factor s 2 in front of the integral in I 3 . 
Finally, write for I 3 




-t\ e -i^^') 2 HV-v) 3 dx „ 

o / ^(x'y+iy—v) 2 

'•t-i IK V<x') 2 -h(jf — Jj) 2 


0 V(0 3 +(2/— v) 2 


•] 


dx" I (62) 


Now, when n\y~y\ is sufficiently large, 

£ ta.»Mv-il>l=0 (ViFTt <63) 

and 

a p®-* e -i*VCr ") 2 +(»-,) s 

\Kx"Y+(y-n) % ~ 

<M» 

Neglecting the terms in equation (64) compared with the 
terms in equation (63), 

(65) 

Comparison of expressions (57), (61), and (65) shows that 
the contribution of It to the total may also be neglected. 
Introducing then expressions (61) and (65) into equation (53), 
the following approximation is obtained 

if I ~ff t 

RS K* 

T Drff] d*h (66) 


Equation (66) contains the fundamental simplification of the 
kernel G in the region B a *. 

Noting that the factor of dX/dij in equation (66) docs not 
depend on £, there is written further 


r 

J*iM Oy * Oyjx 


'it W 


(67) 


Take for (X) rt its value zero immediately in front of the lead- 
ing edge and for (X) Ir its value — ivA which follows from 
equations (35) and (38). Therewith, and with equations - 
(36), equation (67) becomes 

L, (68) 


By introducing equation (68) into equation (66), the follow- 
ing equation is obtained: 

a.* " ) " ! 

(69) 

Equation (69) represents the final result of the present 
section. 


REDUCTION OF THE DOUBLE INTEGRAL OVER THE WAKE 

REGION R „* 


In reducing the double integral over the airfoil region R a * 
to a simpler approximate form, use has been made of the 
fact that the span of R„* is appreciably larger than the chord 
of R a *. Evidently this is not the situation for the wake re- 
gion R u * and thus additional considerations are necessary 
for the integral extended over R„. 

Proceed as follows. Write the last integral on the right of 
equation (50) in the form 

;wrr*-r*> (7o) 


The second of the two inner integrals on the right of equa- 
tion (70) may be treated exactly like the integral over f? a *. 
The first of the integrals has the property, as will be seen, 
that it is dependent on a; in a simple explicit form. 

Taking the second integral first, there results, in analogy 
to equation (56), 


(£ A (Ae*«*) e -**G df] dy « 

[‘■J-T" t + 




w »-’ D ]jr“ 


e~ iri d£dy 


(71) 


The remaining integral with respect to $ in expression (71) 
will be introduced in evaluated form in the final collected 
form of the results. 

Next the remaining integral in equation (70) is transformed 
as follows 

J A (Ae i '*r) e x -k,y—n)d^dn 

= ^ -^r-(Ae i ' x T) £ |* e~ l,<f+A G{ — cr, y — ij) da-J dy 

J A {Ae^Fuiy-n) dy (72) 


It may be noted that the function F u reduces to the 
function F first introduced by Cicala when M= 0, that is, 
when the fluid is incompressible (references 4 and 5). Com- 
bination of equations (72) and (52) results in the following 
form of ■Par, 


e tw (7(— <r,y— y,K)d, 

jz=ii r 

y—y Jo U- 


— ly — ?zt cr ^-iicVff 2 -H' 2 

v 2 +r 2 


(v^+? + “) 




(•-* / f-tf-i! , \ , 


( 73 ) 
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Now equations (71) and (72) are combined in accordance 
with equation (70) in order to obtain the following approxi- 
mation to the double integral over 


]j § ^^ eivXT ) e ~ iri G di-d n ™ 

i£,* 

e ~ ivX f (A + 
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The second integral on the right of equation (75) becomes 

I" [-*? fef «■«*-«+ 


j —’u — t, (W“')^[4 t a»(*-,D+ 


(74) 


THE APPROXIMATE INTEGRAL EQUATION OP THE PROBLEM 

Equations (74) and (69) axe substituted into equation (50) 
and the following approximate integral equation is obtained: 

1 £ x t 

9 (x, y)= \&y)K(x- f; K )d£+ 

iw r* 

— A(y)e ira rj^ e- { ^K(x—^ «)d£4- 

£44 ^ ^ [4-4" Nf • +■ 

W2/-17 1)J di? - ^ e (Ae^Fu d n 

(75) 

Equation (75) represents in preliminary form the result to 
be obtained in this report. This result is reduced to a 
somewhat simpler form as follows. 

Set 

Ae <w! r=Q (76) 

and introduce a new dimensionless span wise coordinate y* 
defined by 

„*= t 


Wl -M 2 


-(77) 


If b represents the semichord at midspan and sb the semi- 
span, then it follows from the definitions of x and y in 
equation (19) that the coordinate y * assumes values in the 
interval (—1,1) only. 

Further, set 


x L (y)=x L *(y*), x T (y)=x T *(y*) 


(78) 


and note, that x L *(0)= — 1, z r *(0) = l. For a rectangular 
lifting surface, —Xj^—Xt*—! throughout. 

From equations (51) and (65), there follows for the first 
integral on the right of equation (75), 

fZ Ul,y) «««*-«+■ 

( 79 ) 


W-t (80) 

The third integral on the right of equation (75) becomes, 
with k from equation (24), 

jli [« T «*«»-- >D+ 

<£* <*8 \y*—y*\ r , iM iir If 

J-1 dv*Y=^ Ll-A^J-. T //o + 

W3G 5 y ’J* 1 )] (81) 

Finally, there is written for the fourth integral on the right 
of equation (75), 

f~(Ae‘”T)F M (y~v)dr,= 

(82) 

where the function F M , as defined by equation (73), may be 
written in the following alternate form: 

«^/;-[/::^( y&? +u 0* + 


/: 


-w p-lM Vr S+r* 


r jd 

r 2 +a: 8 


VrM-'f* 
1 


0 


-rff Jrfr + 


\Vt*-+ * 


i (/r J dtr 


(83) 


Now. ’collect equations (82), (81), (80), and (79) and sub- 
stitute the result in equatiou (75). This gives 

j: 




i?o®( 


(Irl)df J d£+ 
fc c -‘'* 4“ [~ 1 fEf 

/_7“ e ^®(if[)df]^+ 


i/1 — M S J ( Ji— A/s 2 X 


4irS-\/T 

Mct 1 "*-’* 1 ) 


+ 
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ffo (2) (!ri)^rJ- 

. - Fy[- ks -= (g*-i?*)~l} d v * (84) 

Vi -M* \_JT=W y 3 


In equation (84) there are introduced dimensionless coor- 
dinates a:*and £* as follows: 


„* 2x—(x T -\-x L ) 

X = 

X T —XL 


(85) 


so that the interval x L ^x&x T goes over into the interval 
— 1 x* g 1 . Further, 

x T —x L =2b* 

^ (x T -\-x£=x m * 

where b* is the ratio of local semichord to semieliord b at 
midspan. Finally, the following quantities are introduced: 

K*=b*K \ 

v*=b*v > (87) 

a*=u/bV 



Note that, for the rectangular pla n form, 6 *=1, x m *=0 and, 
for the elliptical plan form, 6*= Vl— (y*) 2 , z M *=0. 
Equation (84) then assumes the following form: 


g*(x*, y*) = -3$^ *■*(?, V*) [W - f)l df- + 

~ a *(y*) e-"' 1 -* P r*** K*N [k*(x*~ |*)] df-r 

t) 1 


e -.V(r*+x„.) ri dQ * 


4ts Vi • 


( 88 ) 

Comparison of equations (88) and (84) gives for the ker- 
nels N and Sm the following expressions: 


n ( x)=—~ [M (89) 


Sm(.x)- 


ix ksM |a:l 
= 2 Vl-M 2 * 
ks 


[ffi® (Mist) + fli® (lr i) dij- 




Fu(x) 


(90) 


The task from here on is the following. Equation (88) 
must be solved for X*, in terms of g * and a*. This part of 
the problem is exactly as in the two-dimensional theory. 
The function a* is then to be determined by an integro- 
differential equation which is obtained by expressing a* in 
terms of X* in accordance with the definition of a*. This 
part of the problem is similar to the earlier work on incom- 
pressible flow of references 1, 2, and 6. Finally the solution 
of the integral equation for a* must be used to obtain 
expressions for the pressure ‘ distribution at the airfoil, as 
affected by the three-dimensionality of the flow about a 
wing of finite span. 

The results, as expressed by equations (88), (89), (90), and 
(83), include the special case of two-dimensional flow for 
which da*/dij*=0, and the special case of incompressible 
flow for which ilf=0. They also include essentially known 
results on compressible steady flow for which k=0. 


Massachusetts Institute of Technology, 

Cambridge, Mass., May 18, 1948. 
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